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We show that the low-energy electronic structure of arbitrarily stacked graphene multilayers with 
nearest-neighbor interlayer tunneling consists of chiral pseudospin doublets. Although the number 
of doublets in an A''-layer system depends on the stacking sequence, the pseudospin chirality sum 
is always A^. A*'-layer stacks have N distinct Landau levels a.t E = Q for each spin and valley, and 
quantized Hall conductivity a^y ~ ±(4e'^//i)(7V/2 -I- n) where n is a non-negative integer. 



I. INTRODUCTION 

The recent explosiorP of research on the electronic 
properties of single layer and stacked multilayer graphene 
sheets has be en driven by advan ces in material prepara- 
tion method^^is]^ i^y ^j^g unusuaP*^ electronic proper- 
ties of these materials including unusual quantum Hall 
effect^^, and by hopes that these elegantly tunable sys- 
tems might be useful electronic materials. In this paper, 
we demonstrate an unanticipated low-energy property of 
graphene multilayers, which follows from an interplay be- 
tween interlayer tunneling and the chiral properties of 
low-energy quasiparticles in an isolated graphene sheet. 
Our conclusions apply in the strongest form to models 
with only nearest-neighbor interlayer tunneling, but are 
valid over a broad field range as we explain below. We 
find that the low-energy band structure of any graphene 
multilayer consists of a set of independent pseudospin 
doublets. Within each doublet, the bands are described 
by a pseudospin Hamiltonian of the form 

ldj{k) cx [cos(J(/.fc)r" ± sin(J0fc)T''], (1) 

where r" is a Pauli matrix acting on the doublet pseu- 
dospin, k is an envelope function momentum measured 
from either the K or K' corner of the honeycomb lat- 
tice's Brillouin-zon^, k — |fc|, and is the orientation 
of k. The ± sign in Eq. ([T]) assumes the opposite signs 
in graphene's K and K' valleys. Following the earlier 
work on graphene bilayer^, we refer to J as the chirality 
index of a doublet. In the presence of a perpendicular 
magnetic field B, 7ij(fe) yields the J Landau levels at 
E ^ Q imd E ^ Q levels with \E\ cx B''/'^. Taking the 
twofold spin and valley degeneracies into account, the 
number of independent zero-energy band eigenstates at 
the Dirac point (fc — 0) is therefore 8Ad, where Nd is the 
number of pseudospin doublets. We find that, although 
Nd depends on the details of the stacking sequence. 



= N 



(2) 



in an TV-layer stack. It follows from Eq. ([2| that the 
Hall conductivity of an TV-layer stack has strong integer 
quantum Hall effects with plateau conductivities, 



= ±- 



4e^ 



N 
~2 
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where n is a non-negative integer. 



FIG. 1: (Color online) (a) Energetically favored stacking ar- 
rangements for graphene sheets. The honeycomb lattice of a 
single sheet has two triangular sublattices, labeled by a. and 
/3. Given a starting graphene sheet, the honeycomb lattice 
for the next layer is positioned by displacing either a or /3 
sublattice carbon atoms along a honeycomb edge. This re- 
sults in three distinct two-dimensional sheets, labeled by A, 
B, and C. Representative a and /3 sublattice positions in A, 
B, and C layers are identified in this illustration. It is also 
possible to transform between layer types by rotating by 60° 
about one of the carbon atoms, (b) Each added layer cycles 
around this stacking triangle in either the right-handed or 
the left-handed sense. Reversals of the sense of this rotation 
tend to increase the number of low-energy pseudospin dou- 
blets TV_D. In graphite, the Bernal (AB) stacking corresponds 
to a reversal at every step and orthorhombic (ABC) stacking 
corresponds to no reversals. 



II. STACKING DIAGRAMS AND 
PARTITIONING RULES 



When one graphene layer is placed on another, it is 
energetically favorabl^iSl for the atoms of either a or /3 
sublattices to be displaced along the honeycomb edges, 
as illustrated in Fig. [l] This stacking rule implies the 
three distinct but equivalent projections (labeled A, B, 
and C) of the three-dimensional structure's honeycomb- 
lattice layers onto the x-y plane and 2^~^ distinct TV- 
layer stack sequences. When a B layer is placed on an 
A layer, a C layer on a B layer, or an A layer on a C 
layer, the a sites of the upper layer are above the /3 sites 
of the lower layer and therefore linked by the nearest 
interlayer neighbor 7r-orbital hopping amplitude t\^. For 
the corresponding anticyclic stacking choices (A on B, B 
on C, or C on A), it is the /3 sites of the upper layer and 
the a sites of the lower layer that are linked. 
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FIG. 2: (Color online) Stacking sequences and linkage dia- 
grams for A'^ = 3, 4, 5 layer stacks. The low-energy band and 
the Landau level structure of a graphene stack are readily read 
off these diagrams as explained in the text. Shaded ovals link 
a and /3 nearest interlayer neighbors. 



All distinct iV = 3, = 4, and N — 5 layer stacks are 
illustrated in Fig. |2] in which we have arbitrarily labeled 
the first two layers starting from the bottom as A and 
B. The low-energy band and the Landau level structure 
can be read off these diagrams by partitioning a stack 
using the following rules, which are justified in the fol- 
lowing section, (i) Identify the longest nonoverlapping 
segments within which there are no reversals of stack- 
ing sense. When there is ambiguity in the selection of 
nonoverlapping segments, choose the partitioning which 
incorporates the largest number of layers. Each segment 
(including for interior segments the end layers at which 
reversals take place) defines a J-layer partition of the 
stack and may be associated with a chirality J doublet, 
(ii) Iteratively partition of the remaining segments of the 
stack into smaller J elements, excluding layers contained 
within previously identified partitions, until all layers are 
exhausted. Chirality decompositions which follow from 
these rules are summarized in Table |T1 Note that this 
procedure can result in J = 1 doublets associated with 
separated single layers which remain at the last step in 
the partitioning process. 

In applying these rules, the simplest case is cyclic ABC 
stacking for which there are no stacking sense reversals 
and therefore a single J = N partition. In the opposite 
limit, AB stacking, the stacking sense is reversed in every 
layer and the rules imply A^/2 partitions with J = 2 for 
even N, and when N is odd a remaining J — I parti- 
tion. Between these two limits, a rich variety of qualita- 
tively distinct low-energy behaviors occur. For example, 
in the ABCB stacked tetralayer, ABC is identified as a 
J = 3 doublet and the remaining B layer gives a J = 1 
doublet. The low-energy band structure and the Lan- 
dau level structure of this stack, as illustrated in Fig. [3] 
have two sets of low-energy bands with \E\ oc k, fc^. Lan- 
dau levels with \E\ oc B^/^B^/^ and four zero-energy 
Landau levels per spin and valley. All these properties 
are predicted by the partitioning rules. We have explic- 
itly checked that the rules correctly reproduce the low- 



energy electronic structure for all stacking sequences up 
to = 7. Because each layer is a member of one and only 
one partition, the partitioning rules imply the chirality 
sum r ule in Eq. (|2|. 
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TABLE I: Chirality decomposition for N — 3,4,5,6 layer 
stacks. 
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FIG. 3: (a) Band structure and (b) low-energy Landau levels 
of tetralayer graphene with ABCB stacking. These bands 
were evaluated for nearest intralayer neighbor hopping t = 3 
eV and nearest interlayer neighbor hopping t± = O.lt. 
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III. DEGENERATE STATE PERTURBATION 
THEORY 

We start from the well-known J = 1 massless Dirac 
equatiorfS k ■ p model for isolated sheets, 



T-iMoik) = - 



VTT^ 
VTT 



(4) 



where tt — hv{kx + iky) and v is the quasiparticle ve- 
locity. In the presence of an external magnetic field, tt^ 
and TT are proportional to the Landau level raising and 
lowering operators, so that Eq. Q implies the presence 
of one macroscopically degenerate Landau level at the 
Dirac point for each spin and valle y, and therefore, to 
the iV = 1 quantum Hall effeciP^ of Eq. An N- 

layer stack has a two-dimensional band structure with 
2N atoms per unit cell. The Hamiltonian can be written 
as 



H = Ti.^ + Til 



(5) 



where Ti.± accounts for interlayer tunneling and 7i|| for 
intralayer tunneling. Tiy is the direct product of mass- 
less Dirac model Hamiltonians Hmd for the sublattice 
pseudospin degrees of freedom of each layer. We con- 
struct a low-energy Hamiltonian by first identifying the 
zero-energy eigenstates of H± and then treating Tiy as a 
perturbation. 

Referring to Fig. [2] we see that H± is the direct prod- 
uct of a set of finite-length one-dimensional (ID) tight- 
binding chains and a null matrix with dimension equal 
to the number of isolated sites. An elementary calcula- 
tion gives the following eigenvalues and eigenvectors for 
a chain of length M: 



Er = 2fj^ C0s6'r, 



(6) 



ar = 



M + 1 



(sin Or , sin 26*^ , • • • , sin MOr ) , 



where r = 1, 2, . . . , M is the chain eigenvalue index and 
Or ~ rir/ (M -1-1). Note that odd AI chains have a zero- 
energy eigenstate with an eigenvector that has nonzero 
amplitudes, constant in magnitude and alternating in 
sign, on the sublattice of the chain ends. The set of zero- 
energy eigenstates of 'H± consists of the states localized 
on isolated sites and the single zero-energy eigenstates of 
each odd-length chain. 

The low-energy effective Hamiltonian is evaluated by 
applying leading order degenerate state perturbation the- 
ory to the zero-energy subspace. The matrix element 
of the effective Hamiltonian between degenerate zero- 
energy states r and r' is given bjEl 



(7) 



where n is the smallest positive integer for which the ma- 
trix element is nonzero, and Q = 1 — P, P is a projection 
operator onto the zero-energy subspace. To understand 
the structure of this Hamiltonian, it is helpful to start 
with some simple examples. 



For ABC stacked A''-layer graphene, the zero-energy 
states are the two isolated site states in top and bot- 
tom layers, ai and (3^. The high-energy Hilbert space is 
formed from a set of two-site chains. Because TJy is diag- 
onal in layer index and H± (and hence H^^) can change 
the layer index by one unit, the lowest order at which ai 
and (3n are coupled is n = iV. It follows that 



N 



(8) 



Thus, the effective Hamiltonian of TV-layer graphene with 
ABC stacking has a single J = N doublet^^. 

For AB stacked iV-layer graphene, the high-energy 
Hilbert space consists of a single A'^-site ID chain, exclud- 
ing its zero- energy eigenvalue when N is odd. There is an 
isolated site in each layer which is connected to both its 
neighbors at order n — 2 forming an isolated site chain. 
When N is even, this chain is diagonalized by A^/2, J = 2 
doublets formed between a-sublattice and /3-sublattice 
chain state^I^. When N is odd, the zero-energy chain 
state is mapped to an equal-magnitude oscillating-sign 
linear combination of isolated site states by intralayer 
tunneling at order n = 1, yielding a J = 1 doublet. The 
(TV — l)/2, J — 2 doublets are then formed between a- 
sublattice and /3-sublattice isolated site chain states in 
the orthogonal portion of the isolated state subspace. 

A more complex and more typical example is realized 
by placing a single reversed layer on top of ABC stacked 
N-layer graphene with N > 2. In this case, we obtain 
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(9) 



using a (a^v+i, /^at+i, ai, /3Af) basis. Here, v = VTr/t± and 
is the three-site chain zero-energy state, l3]^^i — 

{Pn+i ~ Pn-i) /V2. The first 2x2 block in Eq. ^ gives 
a J = 1 doublet with a reduced velocity. The J = N dou- 
blet in this instance includes both the {ai,PN) subspace 
contribution and an equal contribution due to perturba- 
tive coupling to the (a^v+i, /S^v+i) subspace. The final 
Hamiltonian is reduced to 



N+1 



N, 



(10) 



where 



-Hi 



-t± 







-t± 



(j.t)W 
(j.)^ 

(11) 

The relationship between the electronic structure of a 
general stack and the partitioning procedure explained 
above can be understood as follows. First, note that a 
partition with chirality J has isolated sites in its terminal 
layers that are coupled at order J in perturbation the- 
ory. In the case of J = 1 partition, the chain opposite 
to the single isolated site always has an odd length and 
provides the zero-energy partner; isolated site to chain 
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coupling therefore always occurs at first order. Next, 
consider the perturbation theory, truncating at succes- 
sively higher orders. When truncated at first order, the 
J = 1 partitions are isolated by higher J blocks within 
which the Hamiltonian vanishes. Each J — 1 partition 
therefore yields a separate massless Dirac equation with 
velocitiefllS that can be smaller than the graphene sheet 
Dirac velocity. When the perturbation theory is trun- 
cated at second order, the Hamiltonian becomes nonzero 
within the J = 2 partitions. The eigenenergies within 
the J — 1 partitions are parametrically larger, and the 
Hamiltonian within the J > 2 partitions is still zero. To 
leading order therefore, the J = 2 partitions are sepa- 
rated, and their isolated states are coupled at the sec- 
ond order in perturbation theory so that each provides 
a J = 2 doublet such as that of an isolated bilayer. If 
two or more J = 2 partitions are adjacent, then their 
Hamiltonians do not separate. In this case, there is a 
chain of second order couplings between isolated states, 
such as those of an even-length AB stack, but the end 
result is still J = 2 doublet for each J = 2 partition. The 
identification between partitions and chiral doublets can 
be established by continuing this consideration up to the 
highest values of J which occur for a particular stack. 
Then, the effective Hamiltonian of any iV-layer graphene 
is as follows: 

Wj^ ^nj,^nj,(g>---^nj,^, (12) 

with the sum rule in Eq. Q. Note that No is half the 
sum of the number of isolated sites and the number of 
odd-length chains. 



IV. DISCUSSION 

The minimal model we have used to derive these re- 
sults is approximately valid in the broad intermediate 
magnetic field B range between ~ 10 and ~ 100 T, 
over which the intralayer hopping energy in field (~ hv/£ 
where i = yJhc/eB ~ 25nm/[B(T)]i/2 is the magnetic 
length) is larger than the distant neighbor interlayer hop- 
ping amplitudes that we have neglected (72 ^ —20 meV), 
but still smaller than t^. For example, if we consider 
ai 03 hopping process in ABA stacked trilayer in Fig. 
[2] the vahd range of magnetic field for the minimal model 
is given by 

h2\<^—^<t^. (13) 

When 72 does not play an important role (in N = 2 
stacks, for example), the lower limit of the validity range 
is parametrically smaller. The minimum field in bilayers 
has been estimated to be ~ by comparing intralayer 
hopping with the 73 ^ 0.3 eV interlayer hopping ampli- 
tude, 

hv^ll < < tj_, (14) 




■0.010 -0.005 0.000 0.005 0.010 ^.010 -0.005 &.000 0.005 fl.010 



FIG. 4: Landau level spectrum at the K valley as a function 
of 73 for an AB stacked bilayer for (a) B = 0.1 T and (b) 
B = 1 T, and as a function of 72 for an ABA stacked trilayer 
for (c) B = 1 T and (d) B = 10 T. Here t = 3 eV, = o'.lt, 
and ujc = eB/mc, with m = t±/2v'^, were used. 

where t;3 — {V3/2)aj:}/h and a is a lattice constant of 
graphene. 

Figure |4] shows the Landau level spectrum at the K 
valley as a function of 73 for an AB stacked bilayer and 
as a function of 72 for an ABA stacked trilayer. In the 
case of the bilayer, the dependence of the Landau levels 
on 73 is weak for B larger than 1 T, whereas in the case 
of the trilayer, the Landau level spectrum still strongly 
depends on 72 for _B = 1 T, but the dependence becomes 
weak for B above 10 T, confirming the above argument. 

In Fig. |5] we plot the noninteracting Hall conduc- 
tivity as a function of the Fermi energy for the four 
distinct four-layer stacks. When electron-electron in- 
teractions are included at an electrostatic mean-field 
(Hartree) level, charge redistribution changes the posi- 
tions of plateaus which are away from the Dirac point, 
but it does not alter the conductivity jumps at the Dirac 
point. If we allow the remote hopping, a small axy = 
plateau appears at the Dirac point for the ABAB staked 
tetralayer due to the energy gap opened by 72. This 
small gap is likely to be important only in very weak 
disorder samples, in which case electron-electron inter- 
actions beyond the Hartree level are also likely to be 
importantP31. The property that the Hall conductivity 
jumps by four units on crossing the Dirac point for ar- 
bitrarily stacked tetralayer graphene is the most obvi- 
ous experimental manifestation of the chirality sum rule 
discussed in this paper. In practice charged multilayers 
{Ep 0) would normally be prepared by placing the 
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FIG. 5: (Color online) Noninteracting Hall conductivity as 
a function of the Fermi energy for all inequivalent four-layer 
graphene stacks when B — 10 T, t = 3 eV, and t± = O.lt. 



many stacks, the A''-fold degeneracy at = persists 
even outside the field range over which the minimal 
model is accurate. Because of the topological charac- 
ter of the quantum Hall effect, we expect the plateaus at 
Wxy\ — (4e^//i)A^/2 to be exceptionally robust in most 
stacks, and that plateaus at the smaller Hall conductivity 
samples will require exceptionally high quality samples. 
We do not anticipate strong quantum Hall effects with 
Wxy \ < {^e^ /h)N/2 unless gaps between the Landau lev- 
els are enhanced by interaction^^ in very high quality'^ 
samples. 

Finally, we note that chiral two-dimensional electron 
system tend toward momentum-space vortex states in 
which charge is spontaneously shifted between layer^^^ 
and that these instabilities are stronger in systems with 
larger J. The present work identifies ABC stacked mul- 
tilayer graphene as the most likely candidate for these 
exotic states. 



system on one side of an electrode and gating. Even 
though gating will redistribute charge and shift energies 
differently in different layers, the Landau level bunching 
we discussed should still be clearly reflected in quantum 
Hall effect measurements. 

In numerical calculations, we have found that for 
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